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The modes of propagation in optical systems which do not possess 
meridional planes of symmetry {nonorthogonal systems) are investigated 
in the case where the effect of apertures and losses can be neglected. The 
fundamental mode of propagation is obtained with the help of a complex 
ray pencil concept. An integral transformation of the field, based on a 
quasi-geometrical optics approximation and a first-order expansion of the 
point characteristic of the optical system, is given; it shows that the complex 
(three-dimensional) wavefront of the fundamental mode is transformed 
according to a generalized "ABCD law." A simple expression is also 
obtained for the phase-shift experienced by the beam. The higher order 
modes of propagation are obtained from a power series expansion of the 
fundamental mode. These higher order modes are expressed, in oblique 
coordinates, as the product of the fundamental solution and finite series of 
Hermite polynomials with real arguments. In the special case of systems 
with rotational symmeti-y, these series reduce to the well-known generalized 
Laguerre polynomials. The theory is applicable to media such as helical 
gas lenses and optical icaveguides suffering from slowly varying deforma- 
tions in three dimensions. Nonorthogonal resonant systems are also in- 
vestigated. An expression for the resonant frequencies, applicable to any 
three-dimensional resonator, is derived. Numerical results are given for the 
resonant frequencies and the resonant field of a twisted path cavity which 
exhibits interesting properties: the usual polarization degeneracy is lifted 
and the intensity pattern of all of the modes possesses a rotational symmetry. 

I. INTRODUCTION 

An optical system, or a resonator, is called "nonorthogonal" when 
it is not possible to define two mutually orthogonal meridional planes 
of symmetry (Ref. 1, p. 240). The helical gas lens 2 - 3 is an example 
of a nonorthogonal lenslike medium. A conventional ring type cavity 
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generally ceases to be orthogonal when its path is twisted, i.e., be- 
comes nonplanar.* 

Let us briefly review the major approaches in the theory of optical 
resonators. The field in a resonator can be expressed exactly in 
terms of known functions only for a few simple boundaiy surfaces. 
No exact solution is available for nonorthogonal systems. However, 
we are interested only in the high frequency operation of large 
resonators. In that limit, the waves have a tendency to follow closed 
curves in the resonator, either clinging to the concave parts of the 
boundary (whispering gallery modes 5 ) or connecting opposite points 
of the boundary (bouncing ball modes). One defines the axial mode 
number as the number of wavelengths existing along such closed 
curves. The nodes of the field in the transverse planes define the 
transverse mode numbers. More insight concerning the mode structure 
and the resonant frequencies can be gained by using a geometrical 
optics approximation, or a paraxial form of the Huygens diffraction 
principle. The geometrical optics approach was developed by Keller 
and Rubinow. 6 It consists of setting up in the resonator a manifold 
of rays tangent to a caustic. The location of the caustic and the 
resonant frequencies are obtained from the condition that the varia- 
tions of the eikonal along three independent closed curves are equal 
to an integral number of wavelengths (or an integer plus one-half 
or one-quarter). This theory, which is analogous to the Born ap- 
proximation of quantum mechanics, gives the exact resonant fre- 
quencies of paraxial modes. The geometrical optics field, when ex- 
tended in the shadow of the caustic by analytic continuation, provides 
an acceptable approximation to the exact field for large transverse 
mode numbers but, for the fundamental mode, it differs vastly from 
the exact field. The caustic line however, does coincide, in two di- 
mensions, with the mode profile. 7 - 8 This geometrical optics method 
has been extended to nonorthogonal resonators incorporating homo- 
geneous media by Popov, 9 who gave an expression for the resonant 
frequencies. Within the paraxial approximation, exact solutions for 
the field can be obtained from the Huygens principle; for that reason, 
the geometrical optics method, in spite of its general interest, will 
not be discussed further in this paper. 

For the case of resonators incorporating inhomogeneous media, 
the Huygens principle must be supplemented by a quasi-geometri- 
cal optics approximation. This approximation consists of assuming 
that a point source at the input plane of the system creates at the 
output plane a field which can be adequately represented by the 
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geometrical optics field. This approximation is generally applicable 
to optical waveguides and resonators if one disregards the effect of 
apertures and assumes that no diffraction gratings or other wave- 
length-dependent scatterers are present. This quasi-geometrical optics 
method provides an integral transformation for the field which is 
equivalent to a partial differential equation of the parabolic type 
(see Section II). The similarity between this parabolic equation and 
the Schroedinger equation has often been pointed out. 10-13 The 
matched modes of propagation in uniform lens-like media with 
hyperbolic secant refractive index laws, for instance, can be found 
in Landau and Lifshits' Quantum Mechanics 1 * [whereas the ray 
trajectories are given in Ref. (1), p. 179]. The more general prob- 
lem of unmatched beams in nonuniform lens-like media corresponds 
to the time-dependent Schroedinger equation with time- va lying po- 
tentials. The adiabatic approximation usually applied to this problem, 
is based on conditions 12 which are too stringent for most optical 
systems. Generalized modes, where allowance is made for a wave- 
front curvature, were introduced by Goubau and Schwering 15 and 
Pierce 16 for the free-space case, in agreement with the theory of 
confocal resonators proposed by Boyd and Gordon. 17 These results 
were extended to orthogonal square law media. 18 ' 19 - 20 The transforma- 
tion of the complex curvature of beams through arbitrary optical 
systems with rotational symmetry and the resonant frequency of 
linear cavities was obtained by Kogelnik. 21 ' 22 Vlasov and Talanov 23 
have observed that, in two dimensions, the phase shift experienced by 
a matched beam in an optical system is equal to the phase of one 
of the two ray-matrix eigenvalues. This result is easily demonstrated 
and generalized to astigmatic orthogonal systems by using a complex- 
ray pencil concept. 4,24 

The generalization to nonorthogonal systems is substantially more 
intricate. Arnaud and Kogelnik 24 have obtained a generalized gaus- 
sian mode of propagation in free space by giving complex values 
to the three parameters which define an astigmatic ray pencil, i.e., 
the position of the focal lines and the angular orientation of one of 
them. This solution can be used to obtain the beam transformation 
in a sequence of thin astigmatic lenses arbitrarily oriented, by 
matching the complex wavefronts at each lens. This method does not 
give, however, a general expression for the phase shift experienced 
by the beam, knowledge of which is essential in studying resonators. 
For that reason, a somewhat different approach is used here, where 
the ray pencil is defined by two of its rays. The field of the funda- 
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mental mode of propagation is obtained (Section III) by allowing 
these two rays to assume complex positions while remaining solu- 
tions of the ray equations. 

The higher order modes of propagation are studied in Section IV. 
They are obtained by application of differential operators related to 
those used in quantum mechanics. An oblique coordinate system 
is introduced which diagonalizes the complex wavefront of the funda- 
mental mode. In this oblique coordinate system, the higher order 
modes can be expressed as the product of the fundamental solution 
and finite series of Hermite polynomials with real arguments. An 
alternative procedure is also given which leads to Hermite polynomials 
in two complex variables. The simple formula for the resonant fre- 
quencies of linear resonators given by Popov ' 20 is shown to be ap- 
plicable to ring type resonators incorporating inhomogeneous media 
(Section V). Finally these general results are applied to a new type 
of optical resonator called "cavity with image rotation" which pre- 
sents interesting resonance and polarization properties (Section VI). 
Numerical results are presented. 

The present theory is limited to paraxial first-order solutions in 
loss-less isotropic media. As indicated before, it is assumed that no 
apertures or diffraction gratings are present in the system, and the 
problem of mode selection is not discussed. The electromagnetic field 
is treated as a scalar quantity and the polarization effects are in- 
troduced only at a later stage; this is permissible within the paraxial 
approximation. Fresnel reflection at surfaces of discontinuity is also 
neglected. 

II. PARABOLIC WAVE EQUATION AND INTEGRAL TRANSFORMATION OF THE 
FIELD 

In this section an approximate form of the scalar Helmholtz equa- 
tion is derived which is applicable to paraxial beams, i.e., to beams 
propagating at small angles with respect to the system axis. It is 
subsequently compared to an integral transformation derived from 
Huygens principle. 

The scalar Helmholtz equation can be written in a x t , x 2 , z rectang- 
ular coordinate system 

dE.dE.dE, 22, \n n /i\ 

—2 + -T-2 + 7X + k n fa , x 2 , z)E = 0, (1) 

dXi dx 2 dz 

where E is a component of the field and n(xi , x* , z) the refractive 
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index of the medium. Let us introduce a reduced field 

^(s, , x a , z) - E(x t , .t, , z) exp jk J n(0, 0, z) dz , (2) 

and neglect the second derivative of \p with respect to z. This approxi- 
mation physically means that only waves propagating in a direction 
close to the z axis are considered. Denoting n(0, 0, z) by n , for brev- 
ity, one obtains 

This equation can be simplified if one introduces the following 
changes of function and variables 24 

* = nU, (4) 

f = ( dz/no . (5) 

•'n 

One obtains 

H + H - 2 >* 5? + k(n - •*• - °- (6) 

Let us further assume that n" — n 2 , is a quadratic form in x t , x a 

n 2 = nl + n n x\ + 2ni 2 z 1 a; 2 + n 22 x\ ■ (7a) 

n,, , n 12 and n 22 are real functions of z since the losses in the medium 
are neglected. The quadratic form given in equation (7a) describes a 
nonorthogonal optical system when the directions of its axes change as z 
varies. In that case, the diagonal term 2n 12 x l x 2 cannot be eliminated by 
rotating the coordinate system about z. We discuss this general case. 
Let us rewrite equation (7a) , for brevity, in matricial form 

n = nl + r-qr (7b) 

where r denotes a column matrix with elements x y , x 2 and 77 denotes a 
2X2 real symmetrical matrix. The sign ~ indicates a transposition. 
Inserting equation (7b) in equation (6), the wave equation assumes 
the form 

£* ^ (v 2 - 2jk y + k 2 f vr y = (8) 

where V J denotes the laplacian operator in the transverse .t, , x 2 plane. 
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It is henceforth assumed that n 2 — n 2 , is small compared with unity. 
Within this (first-order) approximation, the refractive index law, equa- 
tion (7b), becomes 

n ~ w + frjr/(2n ). (9) 

Let us now consider the ray trajectories. A ray (ft is defined at any 
transverse plane z by its position q(z) and by the projection p(z) on 
that plane of a vector directed along the ray, of length equal to the 
refractive index n. q{z) and p(z) are called respectively the position 
vector and the direction vector of the ray. It is convenient to represent 
these vectors by column matrices whose elements are the vector 
components on x v , x 2 . As long as only fixed coordinate systems are 
used, such matrices can be denoted without ambiguity q(z) and p(z), 
or simply q and p. The exact ray equations are (see, for instance, Ref. 1, 
p. 90) 

p = -n VH(r,p), (10a) 

q = n VMr,p), (!0 b ) 

at r = q. In equation (10) the upper dots denote differentiations with 
respect to f, and H(r, p) denotes the Hamiltonian of the system de- 
fined by 

H(r, p) = -{n - pp)*; (10c) 

V denotes the gradient operator in the transverse x v , x 2 plane, and V P 
denotes a gradient operator relative to the p variables. Within the first 
order approximation [equation (9)], equations (10c), (10a) and (10b) 
reduce respectively to 

H(r, p) = -n - (fr,r - pp)/(2n ), (He) 

p = vq, (Ha) 

and 

q = p. (Hb) 

Equations (11a) and (lib) are called the paraxial ray equations. 

Let us now consider two arbitrary rays, (ft and (ft, denned by their 
position and direction vectors q, p and #, p, respectively, and let the 
"product" of these two rays be defined by the scalar expression 

((R; (R) = qp - Jp. (12) 

((ft; (ft) is sometimes called the Lagrange invariant (see Ref. 1, p. 251). 
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It is easy to show that this quantity is independent of f (or z). Indeed, 
applying equations (11a) and (lib) to both <R and (R, and remembering 
that 7] is a symmetric matrix, one obtains 1, 

|:((R;(R) - jj;(0 - b) = & + qp -fr -~® = 0- (13) 

The Lagrange invariant ((R; (R) plays an important role in the present 
theory. Notice that n does not appear explicitly in equations (8), 
(11a) and (lib). It can therefore be assumed, without loss of generality, 
that no as 1. 

The properties of propagating beams are sometimes more easily 
understood by considering the transformation of the field between the 
input plane and the output plane of an optical system described by its 
point characteristic. Let us now choose as optical axis, for generality, 
an arbitrary ray Ct which need not be a straight line nor even a plane 
curve. Let us further define, at a distance z' from an origin 0, a rec- 
tangular coordinate system x[ , x' 2 , whose axes are oriented respectively 
along the principal normal and the binormal to & (see Fig. 1). At any 
given transverse plane, a ray is defined by its position vector q and its 
direction vector p. Let us assume that there is one ray, and only one 
ray which goes from a point r at z = (input plane) to a point r' at 
z = z' (output plane). This assumption implies, in particular, that the 
planes z = and z = z' are not conjugate. The optical length V(r, r') 
of such a ray is called the point characteristic of the optical system. 
As is well known, the direction vectors of a ray can be obtained from V 
by differentiation (Ref. 1, p. 97) 

p = — V1)(r, r'), (14a) 

p> = V'V(r,r'), (14b) 

at r = q, r' = q'. The primes always denote quantities at the output 
plane z = z' '. 

The law of transformation of the field can be obtained from the 
Huygens principle supplemented by a quasi-geometrical optics approxi- 
mation. 28 The Huygens principle states that each point of an incident 
wavefront can be considered as the source of a secondary wave. The 
quasi-geometrical optics approximation consists of assuming that the 
field created at the output plane of the system by a point source at 
the input plane is adequately represented by the geometrical optics 



f Recall also that, for any conformable matrices a and b, (ab)~ = ha and that, for 
any scalar (one element matrix) c, we have £ = c. 
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Fig. 1 — Optical axis of a ring type resonator. T) denotes the point characteristic 
of the system included between two transverse planes, 2 = and z = z'. 

field. These two assumptions allow us to express the field E'(r') at the 
output plane as a function of the field E{r) at the input plane. Within 
the paraxial approximation, we have 



E'(r') = ±\ _1 ft °° E(r)K(r,r')d 2 r, 



where 



(15a) 



(15b) 



K(r, r') = | d°-V/dx { dx' f |* exp [-jkV(r, r')]. 

The term | d 2 V/dXi dx'j |*, where the bars denote a determinant, is 
obtained by recognizing that the power flowing through a small area 
at the output plane is equal to the power flowing in the corresponding 
cone of rays leaving the point source at the input plane, and using 
equation (14a). 

To first order, the quantity S = V — z' is a quadratic form in x x , 
x 2 , x[ , x' 2 which can be written, in matricial notation 

§ = i(r[/r + rVr' + f'Vr + f'Wr'), (16) 

where U and W are 2X2 symmetric real matrices and V is a 2 X 2 
real matrix. Equation (16) can be rewritten, more concisely 



S = Wr'] 



U V 

V w 



■ whm 



(17) 



Introducing equation (16) in equations (14a) and (14b), one obtains 
linear relations between p, p' and q, q' in the form 
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U V 

V w 
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(18) 



It is sometimes convenient to introduce a ray matrix which relates 
q' , p' to q, p. Simple relations exist between [S] and the ray matrix; 
they are given in Appendix A. 

Let us now go back to the integral transformation and observe that, 
if S is a quadratic form [equation (1G)], the determinant 



d*V/dXi dx', | = | d 2 S/dx { dx'; | = | V 



(19) 



is independent of r and r'. This term can consequently be taken out of 
the integral in equation (15). The integral transformation of the re- 
duced field yj/ [$ = E exp (jkz)] becomes 



W) = ±X _I | V |* Jf \[>(r) exp (-jfcS) d 2 r, 

whose kernel is essentially 

K = | V \ h exp (-;feS). 



(20) 



(21) 



Let us show that, in a rectangular coordinate sj^stem, K represents 
the Green function of the parabolic wave equation, equation (8), i.e., 
that 

£ 'K s (v 2 - 2jk ~ + fc^nVjll V \ h exp (-jkS)} = 0. (22) 

The first term in equation (22) can be written, using equation (16) 

V' 2 [| V | J exp (-#§)] 

= | F |* (-j'fcV /2 § - & 2 V'S- V'S) exp (-jkS) 

= \V \ h exp (-i/cS)(-i/cSpui- W - A; 2 V'S- V'S). (23) 

To evaluate the second term in equation (22) one needs to know the 
derivative of S with respect to z' . We have (Ref. 1, p. 97) 



g = -H(r', p') & 1 + (f' v Y - p'p')/2, 



(24) 



where the paraxial approximation of H, equation (lie), has been used. 
Therefore, introducing the expression for p', equation (14b) in equation 
(24), one obtains 



]£ = !?- 1 = (?Vr'- V'S- V'S)/2. 
oz az 



(25) 
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One also needs to know the derivative of V with respect to z'. It is 
obtained by introducing the quadratic form, equation (16), in both 
sides of equation (25) 

. dU , „. dV , , s ,dW . 
r -3-7 r + 2r -77 r + r jt* r 
rfz cfe az 

= f'n'r' - (fV + f'W)(?r + J7r'). (26) 
Equation (26) shows, upon identification, that 

^j=-VW. (27) 

dz 

Therefore (see Ref . 29) 

|r|F|»-*|7|-»^,7|-*|7|**ar(ir*g) 

= -i I 7 I* Spur 17. (28) 

Upon substitution of equations (23), (25) and (28), one finds that 
equation (22) is satisfied. 

Consequently, within the first-order approximation, one may use 
indifferently the parabolic wave equation, equation (8), or the integral 
transformation, equation (20). Most of the demonstrations given in 
the following sections are based on both formulations. 

III. FUNDAMENTAL MODE OF PROPAGATION 

We know that in the high frequency limit, propagating beams 
closely resemble ray pencils. Let us therefore consider first the field 
of such ray pencils, and subsequently see how this solution can be 
generalized to take into account diffraction effects. 

A ray pencil is, in general, astigmatic; it can be defined, in free 
space, as the manifold of rays which intersect two mutually perpen- 
dicular focal lines. At any point, a surface exists, called the wave- 
front, which is perpendicular to all of these rays. The field of ray 
pencils propagating in inhomogeneous media can be written in a 
x l7 x 2 , z rectangular coordinate system 

Efa , x 2 , z) = Ae-' kS , (29) 

where A and S are real functions of Xi , x 2 and z. A is an amplitude 
factor and S is called the eikonal of the geometrical optics field. The 
surfaces S = constant are the equations of the wavefronts associated 
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with the manifold of rays. Let us assume that one of the rays coin- 
cide with the 2-axis and that the refractive index of the medium is 
unity on that axis. Within the first-order approximation, $ = S — z 
is a quadratic form in the transverse variables x x and x-> , whose co- 
efficients are slowly varying functions of z, and A is independent of 
#1 , .T 2 . * can be written, in matrix notation 

*(r, z) m \fn{z)r, (30) 

where n(z) is a 2 X 2 symmetrical matrix which generally depends on z. 
The law of conservation of power dictates that A and n cannot be 
independent; a wavefront with a positive curvature, for instance, corre- 
sponds to a contraction of the ray pencil as z increases, which necessarily 
results in an increased intensity. To express this relation between A 
and n (transport equation) , let us choose any two rays of the ray pencil 
such as (R and (R. Since (R and (R are both perpendicular to the wave- 
front, one has, from equation (30) 

V = VHq) = W> (31a) 

v = v*($) = m, (3ib) 

where V denotes as before the gradient operator in the .r, , x 2 plane. 
Equations (31a and b) can be written more concisely: 



where we have defined 



P = nQ, (32a) 

Q=[qq), (32b) 

P=[pp]. (32c) 

Equations (31) and (32) show that the product of (R and (R, defined in 
equation (12), is equal to zero at any plane 

(<R; (R) = qp - \ V = 0. (33) 

Any ray defined by a linear combination of (R and (R also belongs to the 
ray pencil since its product with either (R or (R is equal to zero. There- 
fore, the one-parameter manifold of rays c(R, e(R + (R, c(R, (R + e(R, 
with < e < 1, defines a tube of rays in the ray pencil whose cross 
section is a parallelogram with sides eq, eq + q, eq and q + e# (see 
Fig. 2). The area of this parallelogram is given by the length of the 
vector product of q and $ 

h = q x q 2 — q 2 q t = \ Q \. (34) 
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Fig. 2 — An astigmatic ray pencil is defined in free space by the manifold of rays 
which intersect two mutually perpendicidar focal lines such as Fi and F2. At any 
transverse plane the intensity of the field is inversely proportional to the square 
root of the area defined by q and q, the position vectors of any two rays of the ray 

pencil ((R and (R). 

Conservation of power requires that A 2 (z)h(z) be a constant. A(z) can 
therefore be obtained from equation (34). Notice that, at a focal line, 
the sign of h(z) changes from positive to negative. Therefore A(z) « 
[h(z)]~ i becomes imaginary. If one insists on keeping A(z) real, a tt/2 
phase shift must be subtracted from S at such points (anomalous 
phase shift). 

The elements of the wavefront matrix n can also be obtained from 
the components of two rays satisfying equation (33). One obtains, 
solving for n equation (32a) 

Mil = (&Pi - q2Pi)h~\ (35a) 

M22 = (91P2 - qiP 2 )h-\ (35b) 

M12 = M21 = (?iPi - QiPi)h~\ 

= (92P2 - q*p2)h~\ (35c) 

The reduced field of the ray pencil is therefore 



4/(r , z; (R, (R) = zkhT* exp ( — 3 ^ *!*) . 



(36) 



where h and /x are given by equations (34) and (35a, b and c) respectively. 
The sign ambiguity in the expression of \f/ can be resolved only by 
counting the number of focal lines along the ray pencil, from some 
reference plane. 
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Let us now show that the field of a ray pencil, as given by equation 
(36), is a solution of the parabolic wave equation, equation (8), i.e., that 



£^(r, z; (R, (R) = ( V" - 2jk — + k'f^r 



h 



I ■ k _ 
5 exp \-]-rfir 



= 0. 



(37) 



The first term on the right side of equation (37) is 
V 2 [rUxp(-j^r>)] 

= /i _i exp f — j^furj X (— jib Spur p — k 2 ?n 2 r). 
The second term is 



(38) 



-'2jk 



h * exp ( -j~rfir) 



dzL h ~ icxp \~ j i ffxr j 

= h~* exp [—j-rnrj X (jkh/h — k 2 rfir). 

Using now equations (34), (32a) and (lib), one notices that 

h/h = Spur n. 



(39) 



(40) 



Differentiating both sides of equations (31a and b) with respect to z 
and using the paraxial ray equations [equations (11a) and (lib)] one ob- 
tains 



(A + n 2 - v)q = 0, 

(m + m 2 - v)q = o. 



(41a) 
(41b) 



Since q and q are generally linearly independent, it results from equa- 
tion (41) that 



A + /* = v- 



(42) 



Upon substitution of equations (38), (39), (40) and (42) in equation 
(37), one finds that the field of a paraxial ray pencil is, as expected, a 
solution of the parabolic wave equation. 

It is important to remark that it has nowhere been specified that 
q, p, q and p are real quantities. The right side of equation (36) therefore 
remains a solution of the wave equation if (R and (R are allowed to be 
complex valued while remaining solutions of the paraxial ray equations. 
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In that c.ise, n(z), whose elements are given by equations (35a, b andc), 
becomes a complex matrix and the exponential term in equation (36) 
describes the intensity pattern of the beam as well as its wavefront. 
As observed before 25 the axes of the constant intensity ellipse do not 
coincide, in general, with the axes of the wavefront surface. It is possible, 
however, to define at any plane an oblique coordinate system in which 
both the real part of n, corresponding to the beam wavefront, and the 
imaginary part of n, corresponding to the beam intensity, are diagonal. 
This coordinate transformation is given at the end of this section and 
used in Section IV to express in a convenient form the higher order 
modes of propagation. 

h(z), given by equation (34), and therefore the amplitude term A (z), 
become complex quantities too. The ±j ambiguity pointed out for the 
case of ray pencils does not exist any more since the phase of A (z) 
changes in a continuous manner along the z axis. 

Let us now consider an optical system described by its point charac- 
teristic matrix [S] and calculate the transformation experienced by an 
incident gaussian beam whose reduced field has the form given in 
equation (36). Introducing this expression in equation (20), one obtains 
a reduced field at the output plane 

W(r') = ±X _1 | V |* /T* 



• jf M exp I -j | [r(U + ju)r + 2? Vr' + f'Wr']j d 2 r. 



(43) 



The integral in equation (43) is easily integrated if one notices that, 
for any nonsingular square matrix m and any comformable column 
matrices r and s one has 

fmr + 2fs = (f + sm~ l )m(r + m~ l s) — tmT l 8. (44) 

Using equation (44) one finds that, if in is a symmetric matrix 

/» A + 00 

// exp [-(fmr + 2rs)] d 2 r = ir \ m | -i exp (SmT's), (45) 

provided the integral is defined, i.e., provided: f (real part of m) r is a 
positive definite form. Substituting 

m = j\(U + n) (46a) 



and 



s = j | Vr' (46b) 
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in equation (45), one obtains a reduced field at the output plane 
*'(r') = ±(-h | U + n I | V IT* 

• exp j-i | f'[W - ?(U + m) _1 V]r'Y (47) 

This field has the same general form as the input field and describes 
a gaussian beam with a wavefront matrix 

n' = W - V(U + nY'V, (48a) 

or, in terms of the ray matrix (see Appendix A) 

M' - (C + D/i)(A + 5m)" 1 . (48b) 

This interesting relation 1, generalizes the "ABCD law" which describes 
the transformation of the complex wavefront in two dimensions. 21,22 
In some applications, it is also of interest to know the phase shift 
experienced by the beam through the optical system. It is given, from 
equation (47), to within II, by the simple expression 

= kz' - i Phase of ( | U + p | | V \~ l ), (49) 

where kz' is the geometrical optics phase shift. Equation (49) reduces 
to the expression given in Ref . 24 in the case of systems with rotational 
symmetry. One also verifies, after a few rearrangements, that the 
amplitude of the beam at the output plane assumes the form given 
in equation (34), i.e., that 

W - -A | U + M | | V I" 1 = rfg - qftl , (50) 

where q' and q' denote the output (complex) ray position vectors. 
Equations (48), (49) and (50) completely define the transformation of 
fundamental gaussian beams propagating along the axis of nonor- 
thogonal optical systems. 

These solutions are easilygenoralized to the case where the axis of 
the incident beam is a ray (R(q, p), distinct from_the system axis. Let 
yf/{r, z) denote the field of an arbitrary beam and (R denote an arbitrary 
ray; one can show 24 that 

*(r, g;<R) = 0>(r - q, z) exp [-jk(rp - fqp)} (51) 

is a solution of the parabolic wave equation, equation (8). An equivalent 

' The transformation of the complex curvature of gaussian beams through non- 
orthogonal systems has been given before (in a very complicated form) by 
Y. Suematsu and H. Fukinuki. 3 " Equation (48b) can alternatively be obtained 
without integration by writing down the laws of transformation of (real) astig- 
matic ray pencils, as suggested in Ref. 25. 
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result can alternatively be obtained from the integral transformation, 
equation (20), by introducing a change of variables. According to 
equation (51), a general form for the propagation of gaussian beams is 
obtained by introducing the expression for the field obtained before 
[equation (36)] into equation (51) 

^(r, z; <R, (R) = /T* exp \-j | [(? - ~q)»(r - q) + 2fp - ~qp)y (52) 

Notice that q and p need not be real for the right side of equation (52) 
to satisfy the parabolic wave equation. It is_ merely required that they 
satisfy the paraxial ray equations. When (R assumes complex values, 
however, it cannot be interpreted any longer as a beam axis. Such 
solutions, with (R complex, are of interest to generate higher order modes 
of propagation, as shown in the next section. 

Let us now show that the fundamental mode of propagation can be 
written in a form resembling the form obtained in the case of orthogonal 
systems. This can be done by introducing, at each transverse plane, a 
coordinate system in which p. is diagonal. 

The reduced eikonal $ can be written 

$ = \ffir = §(fM> + jrp'r), (53) 

where fi and v denote the real and imaginary parts of n, respectively. 
Two quadratic forms such as ffxr and f^r can be simultaneously 
diagonalized if a proper (generally oblique) coordinate system is intro- 
duced. 31 The explicit expression for this transformation is not necessary 
here, because we are interested only in the general form of the field; 
it is given in Appendix C. To deal with oblique coordinates, it is con- 
venient to introduce a tensorial notation. The expression for the scalar 
product of two real vectors* q and p in oblique coordinates assumes the 
form 

q-p = q'Pi i ( 54 ) 

where g 1 and q denote the (contravariant) components of q, obtained 
by drawing lines parallel to the axes from the tip of the vector q as 
shown in Fig. (3), and where Vi , V* denote the (covariant) components 
of p, obtained by drawing lines perpendicular to the axes. For brevity 
the summation sign over repeated indices is omitted. 



tThe following relations are also applicable to complex vectors since such 
vectors can be defined as linear combinations V r + jV t of two real vectors V r 
and V«. 
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INTENSITY = CONSTANT 



Fig. 3 — This figure represents the oblique coordinate system, defined in the X1X2 
transverse plane, which diagonalizes both the real and the imaginary parts of the 
wavefront (represented schematically by ellipses). The contra variant components 
of the position vector q, and the covariant components of the direction cosine 
vector p are also represented. It is assumed that the unit vectors of the coordinate 
system have a unit length. The index in the rectangular coordinate system is placed 
at a lower position only to distinguish it from the oblique coordinate system. 



The reduced eikonal $ is now written 



(55) 



where x',i = 1, 2 (or x', j = 1,2) denote the contravariant components 
of a position vector r, and n denotes a twice covariant tensor. With this 
notation, equations (31) and (32) are valid in any coordinate system. 
Therefore, in the coordinate system in which /x is diagonal (n 12 = 0), 
we have 



Pi = Mil? , 
p 2 = M229 2 , 

Pi = Mil?', 



(56a) 
(56b) 
(56c) 
(56d) 



Let us set 



Mll = 2i = Ei = C, - 2jk~W 



M22 = 2 = ^i = C 2 — 2jk l w' 2 



(57a) 
(57b) 
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In equation (57), C\- and — 2k~ 1 wJ 2 denote the real and imaginary 
parts of na , respectively. The reason for the notation 2/c~ w~ 2 is that 
Wi represent the beam radii along the coordinate axes, as shown in 
equation (59). 

In the new coordinate system (with base vectors of unit lengths), 
the area of the parallelogram constructed on the vectors q and q is 

sin 7 iciq — q 2 q) (58) 

where y is the angle between the two coordinate axes. The field of the 
fundamental mode, equation (36), can consequently be rewritten 

f 00 (x\ x 2 , z; <R, (R) = (sinyYHq'q ~ ffY)"* 

•exp [-[(z 1 M) 2 + (z Z M) 2 ]} 

• exp j-i | [&&? + C 2 (z 2 ) 2 ]}. (59) 

The first exponential term in equation (59) describes the beam in- 
tensity pattern and the second one describes the wavefront of the beam. 
In the special case where the lens-like medium is orthogonal, one 
may choose (R and (ft in two mutually orthogonal planes. Assuming that 
these planes coincide with the x x z and x 2 z planes, respectively, we have 

q 2 = p 2 = qi = A = 0, (60) 

and equation (59) reduces to the known form (see, for example, Ref . 24) 

|~ k "I 

iM^i . #2 , 2; (R, (R) = fff exp — (.-Ci/m'i) 2 — j g C1Z1 

T 2 -fc 2~\ 

•q 2 ' exp —(x 2 /w 2 ) — j-^C 2 x 2 • 



(61) 



IV. HIGHER ORDER MODES OF PROPAGATION 



Two procedures are given in this section to obtain the higher order 
modes of propagation. One is based on the power series expansion of the 
field of off-set gaussian beams and the other is based on the application 
of differential operators on the fundamental mode. These two methods 
can be shown to be equivalent. They lead however to two different 
representations of the field, one in terms of Hermite polynomials in 
two complex variables and the other in terms of finite series of ordinary 
Hermite polynomials. Both representations are of interest. 

It has been shown in the previous section that the field of a gaussian 
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beam propagating along the axis of an optical system is fully described 
by two complex rays, denoted (R and (R, satisfying equation (33). This 
solution of the wave equation ^an be generalized to include the case 
where the beam axis is a ray (R. It was pointed out also that (R may 
assume complex values provided its position and direction vectors 
(q, p) remain solutions of the ray equations. Let us define (R as a linear 
combination of the rays (R* and (R*, conjugate of (R and (R, respectively. 
We have 



q = a x q* + oc 2 q* = Q*a, 
p = a,p* + a 2 p* = P*a, 



(62) 
(63) 



where a x and a 2 are two arbitrary parameters. Introducing these ex- 
pressions in equation (52) one obtains 

iKr, z; (R, (R; a, , a 2 ) = ^ 00 (r, z; (R, (R) X exp (ay - \ava), (64) 

where ^ 00 (r, z; (R, (R) denotes the fundamental mode field and where we 
have defined 



Lota. 



v n P 12 

.''12 "22J 



= -abflfVO*, 



// ■ 



L2/ 2 J 



= -2kQ*n'r. 



(65a) 



(65b) 



(65c) 



Notice that i» is a symmetric matrix as a result of equation (33). 

One now observes that the exponential term in equation (64) is the 
generating function for Hermite polynomials in two variables 32 



exp (ay — hava) = 2 "^ ~f H mn (v l y,p), 



„=o mi ni 
where the polynomials II mn liave the form 

1 m(m - 1) 



H mn (yi , y 2 ; v) = y"yl - 



1 



"ll2/l «/2 



. mn ,„_,„_, . 1 n(n — 1) „, n _ 2 ~| 

+ "J" ".2//, //2 + g j " "22?/,Z/2 J + ' ' 



(66) 



(67) 
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and, for m + n ^ 3 

#00 = 1> 

#10 = Vl , 

#oi = y 2 , 

#20 = yl — "ii i 

#11 = 2/l2/2 — "l2 ) /gg\ 

#02 — Vi ~ "22 , 
#30 = V\ — %V\\Vl » 

#21 = 2/?i/2 — 2b» H yj — v tl y 2 , 
#12 = VilA — &utih - "222/1 , 

#03 = yl — 3^222/2 • 

Each coefficient in the expansion of yp(r, z; (R, (R; ^ , a 2 ) in power 
series of a, , a 2 is necessarily a solution of the wave equation since a x 
and a 2 are arbitrary numbers. New solutions of the wave equation are 
therefore obtained in the form 

^ mn (r, z; (Jt, (R) = tfoo(r, z; <R, ^)# mn (Q* _1 r; *)• (69) 

It is demonstrated in Appendix D that this set of solutions forms an 
orthogonal system, provided the condition ((R; (R*) = is satisfied 
[in addition to equation (33)]. The fact that y t and y 2 are complex does 
not raise any particular difficulty in calculating H mn (y x , y 2 ; v) from 
equations (67) and (68). This prevents us, however, from identifying 
y l and y 2 with real coordinates. It is important to notice that multiplica- 
tion of (R by a factor X (i.e., q —» Xg, p — > Xp) and (R by a factor 
X (q — > %q, p — > Xp) leaves essentially unchanged the field given in 
equation (69); it is merely multiplied by a constant. This property 
results from equation (65) and the general form of # mn given in equa- 
tion (67). Consequently (R and (R need to be defined only to within 
constant factors. 

In the special case where the optical system is orthogonal, one may 
choose (R and (R in two mutually perpendicular meridional planes co- 
incident with the x v z and x 2 z planes respectively {q 2 = p 2 = 0, 
£j = p, = 0). The matrix v becomes diagonal and the Hermite poly- 
nomials in two variables reduce to a product of two Hermite poly- 
nomials in one variable. To within a constant one has, in that special 
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case 



H mn (y x , y 2 ;v) = ^(^H^y^B^y^ (70) 

where H k {x) denotes a Hermite polynomial in one variable of order k 
(as defined in Ref. 33). Using equation (57), the right side of equation 
(70) can be written, to within a constant 

(q\/qr /2 m/^) n/2 H m (2 i x 1 /w 1 )H n (2 i x 2 /w 2 ) > (71) 

in agreement with previous results. 11 

The procedure just described for obtaining new solutions of the 
wave equation can be applied to an arbitrary field ^(r, z) . The coefficients 
of the power series expansion are obtained in that case by repeated 
differentiation. If one calculates the coefficients for the few first orders, 
one finds that they assume the form 

* mn {r, z; &, &) - A m ((R*)A n ((R*)*(r, z), (72) 

where A((R) and A((R) are differential operators defined by 

A((R) = pr - jk~ l qV, (73a) 

A((R) = ~p- - jfc-^V. (73b) 

It is not difficult to show, using equation (33), that these two operators 
commute with one another. For generality, let us demonstrate equa- 
tion (72) on the basis of the integral transformation, equation (20). f 
Let 4f(r) and V(r') denote fields at the input and output planes, 
respectively, of an optical system described by its reduced point charac- 
teristic S. Let us prove that a field A((R)^(r) is transformed into 
A'((R')^'(?'') at the output plane, i.e., that 

(pY - Jk-'q'V'^fJ +(r) exp (-jfcS) dVj 

= // +C ° [(Pr ~ jk-'qVHir)] exp (-jkS) d 2 r. (74) 

Notice that the constant term ±X _1 | V |* in equation (20) can be 
dropped. The primes in equation (74) refer as before to quantities taken 
at the output plane. 

Using equation (16), one finds that 

V exp (-jkS) = -jk(?r + Wr') exp (-jfcS). (75) 



t Alternatively one can show that the operator A((R) [or A((R)] commutes with the 
wave equation operator, equation (8). This result has been obtained before by 
Popov 26 for a special form of A((R). 



2332 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1070 

Therefore the left side of equation (74) can be written 

// " \P'r' - q'{Vr + Wr')]^(f) exp (-jfeS) d 2 r. (76) 

To evaluate the right side of equation (74), notice that, for any function 
F(Xi , x a ) which tends exponentially to zero as i, , x 2 — * ± ■» , one has 

[J VFfo , x 2 ) dx t dx a = 0. (77) 

Therefore, setting 

F{x, , x 2 ) = *(r) X exp (-jkS) (78) 

in equation (77), one obtains 

[J VHr) X exp (-;fcS) dS 

= -ff (-;AVS)^(r) X exp (-;fcS) d 2 r. (79) 

Using again equation (16) to evaluate VS, the right side in equation (74) 
becomes 



ff" [W + q(Ur + Vr')]$(r) exp (-j/cS) d 2 r. 



(80) 



The identity of the two terms in brackets in equations (76) and (80) 
results from the ray equations, equation (18). 

The property established for A((R) clearly holds true also for the 
operator A m ((R) corresponding to m applications of A((R), and for the 
operator A n ((R) associated with another ray (R. 

When applied to a gaussian beam (denned by 01, (R), the operators 
A (01) and A((R) give a result identically equal to zero. Higher order 
modes are obtained, however, if one considers the operators associated 
with the conjugate rays (R*, (R*. One therefore calculates 

^ mn (r, z; Gt, (R) = A m ((R*)A n ((R*)iM?', *J &> &)• (81) 

To give a convenient form to the right side of equation (81), let us 
write down explicitely in tensorial notation (see Section III) the operator 
A((R*), denned by equation (73 a) 

A((R*) = p*iX* - j/r y*v, 

- (pW - >ry* |t) + fax* - jk-y* £,). (82) 

Using equations (82) and (59) and the relation 33 
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j- x H k (x) = 2xH k (x) - H k+l (x), (83) 

where H k (x) denotes a Hermite polynomial of order k, one finds that 
m*){H^2 i x l /w 1 )H n (2 i xVw 2 )Mx\x\z;(Sl,a)} 

= fao(x\ x 2 , z; (R, (R)[g 1 *Mff- +1 (2VM)Hn(2V/w 2 ) 

+ g 2 */«'A(2V/w 1 )ff, +1 (2VM)]. (84) 

A similar relation holds for A((R*). These two relations show, by recur- 
rence, that the field of the mode m, n can be written 

if/ mn (x\ x 2 , z; (R, (R) - A m (CR*)A"(m*)^oo(a: 1 , x 2 , z; (R, (R) 

= M* 1 !**! *;«i«) X [q 1 */w l H(2 i x l /w l ) 

+ q 2 */w i H(2*x 2 /w 2 )T X [f*/w l H(2 i x x /w 1 ) 

+ ? 2 *Mtf(2*.r 2 M)] n , (85) 

where the convention is made that, a/fe? - multiplication of the two 
binomials, H k (x) actually represents a Hermite polynomial of order 
k: H k (x). This form of the field shows that the higher order modes of 
propagation can be obtained by multiplying the fundamental solution 
by a finite series of Hermite polynomials in one real variable. 1 " Since 
q /q and q/q are generally complex, the wavefronts are different for 
each mode. It is shown in the next section that q 2 /q l and q/q happen 
to be real, however, at the end mirrors of linear resonators. From this 
observation, it results that the wavefronts of all of the resonating 
modes generally coincide with the end mirror surfaces. 

Another special case of interest is the case where q 2 /q and q/q are 
both equal to j. This happens in the case of systems with rotational 
symmetry, such as the "cavities with image rotation" which are investi- 
gated in Section VI. 

V. NONORTHOGONAL RESONATORS 

We are concerned in this section with the resonant fields and the 
resonant frequencies of nonorthogonal resonators. Ring-type resonators 



t In the case where m = (or n = 0) it is not difficult to show that the two 
expressions given for the mode ma in equations (69) and (85), respectively, 
coincide^ This result can be obtained by writing equation (69) in the coordinate 
system in which y. (but not necessarily v) is diagonal and using an expansion 
formula [equation (22), p. 371 of Ref. 32] for H mn and a condensation formula 
[equation (31), p. 345 of Ref. 321 for the right side of equation (85). In the 
general case a direct comparison of the two expressions appears to be difficult. 
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being conceptually simpler than linear resonators, their properties are 
considered first. A ring-type resonator is essentially a section of wave- 
guide closed on itself. An optical beam is a mode of the resonator, if, 
after a round trip, its field reproduces itself exactly. 

The general form of the solutions obtained in the previous sections 
(Sections III and IV) is preserved as the beams propagate through 
an optical system. In general, however, the field distribution at the 
output plane of a section of waveguide does not coincide with the 
field distribution at the input plane [see, for instance, the transforma- 
tion law, equation (48b) for the fundamental mode]. By a proper 
choice of the mode parameters it is possible, however, to achieve 
coincidence between the fields at the two planes (except, perhaps, 
for a constant phase factor). In that case, the beam is said to be 
matched to the section of waveguide considered. Clearly, such a beam 
would also be matched to a sequence of identical sections, forming 
a periodic waveguide. For the fundamental mode, the matching con- 
dition can be obtained by specifying that // = /x in equation (48b) 
and solving for /*.. However it is more convenient to look first for 
rays which reproduce themselves after a round trip in the system 
(except for a constant factor) and calculate the wavefront matrix /x 
associated with these rays. Such rays are called eigenrays; they are 
always complex in the case of stable resonators. 

To obtain the eigenrays, let us replace o/ and p' by \q and 
Ap, respectively, in equation (18). One obtains the relations 

- v = (U + \V)q, (86a) 

p = (A" 1 ? + W)q, (86b) 

and, by addition and subtraction 

= (U + W + XV + \~ l V)q, (87a) 

p = ±(W - U + A -1 ? - \V)q. (87b) 

Equation (87a) actually represents a system of two homogeneous 
linear equations which admit a solution only if 

| U + W + W + X _1 F | = 0, (88) 

where the bars denote a determinant. Equation (88) can be rewritten 
as a second-degree equation in (A + A -1 ) as shown previously 27 for a 
special case. One obtains 

| V | (A + A" 1 ) 2 + [V 11 K M + K U V 22 - K i2 (V 12 + 7 21 )](A + A" 1 ) 

+ \K\- (V 12 - 7 21 ) 2 = 0, (89) 
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where we have defined 

K=U + W. (90) 

The resonator is stable when the solutions of equation (89) for 

(X + A"')/2 = cos 6 (91) 

are real and are in the range —1 to +1; this is assumed henceforth. 
In that case, two real characteristic angles, denoted d and 0, are ob- 
tained, the two other characteristic angles being clearly — and — 0. 

If one introduces one of the four eigenvalues X = exp (jd), X = 
exp (jd), X* = exp (—jd) or X* = exp (—jd) in equations (87a) and 
(87b), one obtains (to within arbitrary constants) the components of the 
four eigenrays denoted respectively (R, (R, (R* and (R*. Let us show that 
the product of (R and <R [defined in equation (12)] is equal to zero. 

Since ((R; (R) is invariant, one may choose a reference plane along the 
path where the matrix V is symmetric. At such a plane, equations (87a) 
and (87b) assume the form 

- [U + W + (X + X" 1 )^, (92a) 

p-l\W-U+ (X' 1 - \)V]q. (92b) 

Since both U + W and V are symmetric, one has 31 

qVq = ~qVq = 0, (93) 

provided the absolute values of 6 and 6 are distinct. Therefore 

((R;(R) = qp - ~qp 

= \q[W - U + (X- - X)7]$ - \q\W - U + (X" 1 - \)V]q 

= MX"' - %)qVq- - i(X~' - \)~qVq = 0. (94a) 

One also has, replacing X by X -1 and/or X by X -1 

((R*;(R*) = 0; ((R;(R*) = 0; ((R*;(R) = 0. (94b) 

Therefore, according to the results of Section III, each pair of eigenrays 
in equation (94) defines a gaussian beam. The choice between the four 
pairs of eigenrays can be made by giving either a positive or a negative 
sign to 6 and 0. It is made in such a way that the imaginary part of y. 
is a negative definite form. This ensures that the pow r er carried by the 
beam in finite. After traversing a period of the optical system, the 
position q and the direction p of (R become q exp (jd) and p exp (jd) 
respectively. Similarly, q and p become q exp (jd) and p exp (jd) . Equa- 
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tions (35a through c) and (34) show that m assumes its original value 
after a period (round trip) in the optical system; h, however, is multi- 
plied by exp [j{6 +0)]. The field of the fundamental gaussian beam 
denned by (R and (R consequently reproduces itself after a period except 
for an additional phase shift equal to kL — (0 + 0)/2, where L denotes 
the period (round trip) path length. 

To clarify the above discussion let us observe that the modal matrix 

«* «1, (94c) 

\JkP* P_ 

where Q and P were defined before in equations (32b and c), is itself a 
ray matrix, i.e., satisfies equation (112). As shown in Appendix D, the 
imaginary part of \i can be written -(k~ 1 /2)(QQ*)~ 1 ; this is clearly a 
negative definite form, as required. It can also be shown, using equations 
(111), (16) and (21), that the mode generating function \ft(r; a) given in 
equation (64) is precisely the output field created by a point source 
located at the input plane of a (lossy) optical system whose ray matrix 
is the modal matrix, equation (94c). 

Considering now the form, equation (85), obtained for the higher 
order modes, it appears that the operators A m and A n are responsible for 
an increase of the phase of $ mn equal to — md — nd. Therefore, the 
general expression for the resonant frequencies is 

k tmn L = (m + |)0 + (n + 1)0 + 2ftr, (95) 

where I is an integer defining the number of wavelengths along the 
system axis. This result was obtained by Popov 9,20 for the special case 
of linear nonorthogonal resonators incorporating homogeneous internal 
media. It is shown here to be applicable to the general case. 

Let us now investigate the case of linear cavities (cavities with 
folded optical axis). It is convenient to replace the two curved end 
mirrors of such resonators by plane mirrors and thin lenses, and take the 
reference plane at one of the end mirrors. In a round trip along the 
folded optical axis two optical systems are encountered which are 
mirror images of one another. It is shown in Appendix B that the point 
characteristic matrix assumes in that case the simple form 



[S] = 



(96) 



U V 

y u_ 

where both U and V are real and symmetric. The characteristic equa- 
tions (92a) and (92b) become simply 

(U + cos 9V)q - 0, (97) 
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p = —j sin 9Vq, 



(98) 



where 6 is the characteristic angle. 

Let q, p and q, p denote two solutions of equations (97) and (98) 
(eigenrays). Equation (97) shows that the ratio of the two components 
of q and q, q 2 /q\ and q 2 /qi , respectively, are real (in any coordinate 
system) since, for a stable resonator, the solutions and of the charac- 
teristic equation 



u + cos ev I = o 



(99) 



are real. One also observes that the wavefront matrix n is imaginary. 
This result shows that, at the end mirrors of linear resonators, the wave- 
front of all of the modes coincide with the mirror surfaces, except 
perhaps in some cases of degeneracy. Since U and V are symmetrical, 
one further notices that 31 



qVq = 0, 



(100) 



provided the absolute values of and are distinct. Therefore, from 
equation (98), 



qp = fiq* = 0. 
This relation is useful in checking numerical calculations. 



(101) 



VI. CAVITIES WITH IMAGE ROTATION 

As an example of application of the general theory discussed in the 
previous sections, let us calculate the resonant frequencies and the 
resonant field of a new type of optical cavity that one may call "cavity 
with image rotation." 

Consider a nonplanar closed path (see Fig. 4) and let ft be the rotation 
experienced after a round trip by rays parallel to the optical axis. 
(The value of ft for a given orientation of the mirrors can be found in 
Ref. 4.) The case where the optical system has a rotational symmetry is 
of particular interest. Let [° b d ] be the 2X2 ray matrix of the optical 
system with rotational symmetry introduced along the path. The round 
trip point characteristic matrix of the resonator is, in rectangular 
coordinates 



[S] = b- 



a 





— cos ft 


-sin ft 





a 


sin ft 


— cos ft 


— cos ft 


sin ft 


d 





—sin ft 


— cos ft 





d 



(102) 
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(PATH LENGTH L 



, IMAGE ROTATION Si 



Fig, 4 — A cavity with image rotation is represented. It incorporates a lens and 
four plane mirrors which define a nonplanar path. As a result of the twist of the 
path, this resonator is nonorthogonal. When the lens is astigmatic, the resonating 
modes do not exhibit the same patterns as in the case of more conventional 
cavities. 

Equations (89) and (102) show that the characteristic angles are simply 

d = d + 0, (103a) 

6 = O - a, (103b) 

where we have defined 

cos O = (a + d)/2. (104) 

The resonant frequencies are therefore given, from the general rela- 
tion equation (95), by 

k tmn L = (to + n + l)0o + (» - n ± 1)0 + 2(t. (105) 

The additional term ±Q in equation (105) is to be introduced when 
polarization effects are taken into account. It has been assumed that 
the mirrors are perfect conductors, even in number, and that the medium 
is isotropic. In that case the polarization vector experiences the same 
transformation as an image, 4 i.e. a rotation Q. The -f- and — signs in 
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equation (105) refer to the clockwise and counterclockwise polarization 
states, whose degeneracy is therefore lifted. 

The eigenvectors q, p and #, p have respectively clockwise and counter- 
clockwise circular polarizations too, as one expects from the rotational 
symmetry of the system; they are independent of the image rotation £2. 
The components of (R(g, p) and &(<?, p) are respectively, to within 
arbitrary constants 

J«(JM) «(«»■_-» (106) 

lp(— sin O , jam 6 ) lp(— sin 6 , — j'sin d ). 

Setting for brevity, 2 i x 1 /w = x, and 2 i x 2 /w = x 2 , where w is the beam 
radius, the mode ^ mn assumes in rectangular coordinates, from equa- 
tion (85), the form 

+ mn = [H{ Xl ) + jH{x 2 )) m [H{x l ) - jH(x 2 )) n ho ■ (107a) 

This expression, being independent of ft, should coincide with known 
forms (see Ref. 15 or 11) which can be written 

(-l) n n!2 m+n Z m - n L;- n (ZZ*)^oo if m ^ n, (107b) 

(-l) m m\2 m+n Z* n - m L n m - m (ZZ*)f 00 if m^n, (107 c) 

where L' p denotes a generalized Laguerre polynomial, and 

Z = Xi+ jx 2 . (108) 

A relation between Hermite polynomials and generalized Laguerre 
polynomials was given before, in a different form, by J. R. Pierce and 
S. P. Morgan (private communication). The identity of the right side 
of equation (107a) and equations (107b) and (107c) is easily demon- 
strated for the special cases n (or m) = 0, and m = n, using well- 
known formulas,* and verified for the first values of m, n. The 
field consequently assumes the same form as in ordinary cavities. A 
rotation CI about the z axis of the beam pattern can be expressed by 
a multiplication of Z by exp (jCl) and consequently, from equation 
(107), by a phase shift (m — n)n, in agreement with equation (105). 
The distinctive feature of cavities with image rotation compared with 
ordinary cavities, in addition to the polarization properties mentioned 
before, is that the intensity pattern of the resonant field has neces- 
sarily a circular symmetry. 

When the optical system introduced along the nonplanar path is 



*Sce Rcf. 33. Notice that a factor 2 2 " is missing in equation (32), p, 195, of this 
book. 
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astigmatic, one must use the general expressions given in the previous 
sections. This case has been studied numerically, using equation (85), 
for the case of a resonator incorporating a single spherical mirror of 
radius R = 6m, operating at an incidence angle of 30° and an odd 
number of plane mirrors. The spherical mirror is equivalent to an 
astigmatic lens of focal lengths /i = 2.6m and / 2 = 3.47m. Assuming 
a round trip path length L = lm and an image rotation Q = 20°, 
one obtains for the point characteristic matrices, from equation (115), 
with d = 1, v = 



U = 



V = 



W = 



0.615 




0.712_ 




-0.94 -0.34] 


0.34 -0.94_ 



The characteristic angles are 



1 
1 



= -13°3, 
= -54°, 



from which the resonant frequencies can be obtained. The components 
of the eigenrays (R and <R are respectively, in a rectangular coordinate 
system 



m 



?(i, -;1.35) 

p(0.19 - jO.66, -0.62 - jO.19), 



~ [0(1, ;0.91) 
lp(0.19 - J0.57, 0.49 + i0.13). 

These two eigenrays fulfill, as expected, the condition qp = gp; they 
define a wavefront matrix 



M = 



0.096 - i0.305 
0.0206 



0.0206 
0.072 - J0.246J 



whose imaginary part is a negative definite form, as required. The 
intensity pattern for the mode \p 2 o is shown in Fig. 5. It is inter- 
mediate between the circularly symmetric patterns observed when 
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RELATIVE 
INTENSITY = 1 




Fig. 5 — This figure represents the constant intensity curves of the TEMao mode 
in a nonorthogonal cavity incorporating a 6m radius mirror with an incidence 
angle of 30° and an image rotation of 20°. The optical axis path length is lm, and 
the wavelength is Ifixn. 

h = h, and the usual orthogonal patterns observed for fi = (see, 
for instance, the TEAL-n mode in Fig. 7 of Ref. 11) . 



VII. CONCLUSION 

It has been shown that, within the first order approximation, the 
solutions of the scalar wave equation can be expressed in terms of the 
solutions of the (simpler) ray equations. The fundamental mode of 
propagation in nonorthogonal media was obtained by generalizing the 
expression for the field of astigmatic ray-pencils. An oblique coordinate 
system has been introduced which reduces this solution to the form 
assumed by ordinary gaussian beams. The higher order modes of prop- 
agation were also obtained; they can be expressed as the product of 
the fundamental solution and Hermite polynomials in one real variable. 
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The results of Popov - 26 for the resonant frequencies of nonorthogonal 
resonators were extended to resonators incorporating arbitrary lens-like 
media and were applied to a new type of cavity which exhibits interest- 
ing resonance and polarization properties. This theory may also be 
useful for special optical waveguides such as the helical gas lenses, and 
for analysis of optical systems which are nominally orthogonal, but 
which suffer from small distortions in three dimensions. 



APPENDIX A 

Relations Between the Point Characteristic Matrix and the Ray Matrix 

It has been shown in the main text [equation (18)] that the direc- 
tion vectors p, p' of a ray at the input and output plane of an optical 
system are related to the position vectors q, of by the following matri- 
cial relation 



-p 

LP' J 



U 



v 



W_ \_q'_ 



= [S] 



-<?'- 



(109) 



where U and W are 2 X 2 real symmetric matrices, V is a 2 X 2 real 
matrix. [S] is a 4 X 4 symmetric matrix which has been called the point 
characteristic matrix. One also sometimes defines a ray matrix, [311] 
which relates the position and direction vectors of a ray at the output 
plane to the values assumed at the input plane 



(110) 



where A, B, C, D are 2 X 2 real matrices. Since, from equation (109), 
only 10 numbers suffice to define the optical system, the elements of the 
4 X 4 ray matrix [311] must be related by 16 — 10 = 6 relations. To 
obtain these relations, let us compare equations (109) and (110). One 
obtains readily 



IV" 


= 


[a b~\ 
c d\ 


[q ] 


^ [an] 


\q\ 


-p'- 




-PJ 




-V- 



U = B~ l A, 

V = -B~\ 

V = C - DB-'A, 
W = DB'\ 

Since U and W are symmetrical one has 

AB - BA = 0, 



(Ilia) 
(Hlb) 
(Ulc) 

(Hid) 



(112a) 
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BD - DB = 0, (112b) 

and, by comparing the expressions obtained for V and P, equations 
(111b) and (111c), and using equation (112b), one finds that 

DA - BC = 1. (112c) 

Equations (112a, b and c) are equivalent to those given by Luneburg. 1 
They effectively correspond to six independent relations. The relations 
inverse of equations (Ilia through d) are 



A = -V- X U, 
B = -V' 1 , 

c = ? - wv l u, 

D = -WV~\ 



(113a) 
(113b) 
(113c) 
(113d) 



APPENDIX B 

Point Characteristic Matrix of a Sequence of Thin Lenses and Mirrors- 
Symmetrical Systems 

The point characteristic matrix [S] of a sequence of thin astigmatic 
lenses and plane mirrors, arbitrarily oriented in space, can be obtained 
in closed form. 

Let us first consider a thin astigmatic lens oriented at an angle v 
with respect to the x x axis of a XiX 2 z rectangular coordinate system, 
with focal lengths /i , / 2 • This lens is followed by a section of free space 
of length d. For generality, one further assumes that the output co- 
ordinate system is rotated by an angle £2 about the z axis. This rotation 
has to be introduced in the case of non planar paths. 4,34 Using the 
expression for the optical thickness of a lens, and the paraxial approxi- 
mation of the length of tilted ra3 r s in free space, one obtains 

V~) 



[s] s 



(H4) 



with 



U = 



1 / cos v sin v 
d \ /, + U 

■ ( l 1 
cosy sin j>l— — - 

Mi h 



W 



cos v sin v I — 

VI U 



cos 2 v sin 2 v 



h 
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V = -d' 



W = (f 



-1 


cos 12 


sin i2 




_— sin £2 


cos 12 


[l 


o| 




„0 


U 





(115) 



These expressions are applicable to curved mirrors under oblique inci- 
dence with little modification since a curved mirror is equivalent to a 
plane mirror and a lens, in the most general case. 34 It remains to calcu- 
late the point characteristic matrix of a sequence of optical systems 
such as the one described by equations (114) and (115). 

The point characteristic matrix [S«] of a sequence of two optical 
systems whose point characteristic matrices are respectively [Si] and 
[S 2 ] is obtained by using equation (18) of the main text, and specifying 
that the rays are continuous at the junction between the two systems. 
One obtains 



[*.] = 



tfi - v l (w l + uj*Vi -v 1 (w 1 + vjr l v a 
L -T,ar, + uf*Vi w 2 - ? 2 <yr l + u 2 y l v 2 \ 



. (H6) 



In the special case where the second optical system is the mirror 
image of the first system with respect to their common plane, [S t ] 
reduces to 



[«.]- 



U, V t 



u - \vw~ l v -\vw~ l f 

-hVW~ l V U - WW 1 ? 



(117) 



where the index 1 has been omitted. Equation (117) shows that, in 
a symmetric system, U, is equal to W t and V t is a symmetric matrix. 
Repeated applications of equation (116) and equations (114) and 
(115) give the point characteristic matrix of an arbitrary sequence of 
lenses or mirrors. 



appendix c 

Diagonalization of a Complex Wavefront 

The need for introducing an oblique coordinate system at each 
transverse plane has been outlined in the main text. Detailed transforma- 
tion formulas are given in this appendix. 

Let e x , e 2 be the base vectors, of unit length, of the original rectangular 
coordinate system, and d , e 2 the base vectors, also of unit length, of a 
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new coordinate system. f The e, , i — 1, 2 are linearly related to the 
e, ,j= 1, 2 by 

e, = 5(e y , (118) 

where 5j is the mixed tensor which expresses the coordinate transforma- 
tion. The reduced eikonal <J> is a complex quadratic form which was 
written in the original coordinate system [equation (53)] 

<*> = if/xr = \r(n r + i/x*>, (119) 

where n' and ju' are real symmetric matrices. 

By stipulating that, in the new coordinate system, the off-diagonal 
terms of // and p* are both equal to zero, one obtains the transfor- 
mation [8] which diagonalizes $ 

(1 + vT 1 ' v(l + vT h ~ 



[«J - 



*1 5? 



(120) 



where 



u - (c/a)i> = [-b + (b 2 - 4ac)*]/2o, (121) 

a = M11M12 — M11M12 , (122a) 

6 = mIiM22 — m!iM22 , (122b) 

c s ii22Mi2 — M22M12 • (122c) 

The law of transformation of the contravariant components of a vector q, 
denoted respectively q in the old system and q* in the new system is 
(omitting the summation sign) 

f/ = SJq''. , (123a) 

This relation is also applicable to the coordinate x' 

x' = 8)x'. ,. (123b) 

The covariant components of a vector p, denoted 7?, in the old system 
and p.- in the new system, transform according to the inverse relation 

P. = «JP, • (124) 

Expressions for the new components of y. are derived in the main text. 



t Quantities relative to the new system are denoted by bold face" letters in this 
Appendix. Ordinary letters are used in the main text, where there is no risk 
of confusion. 
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APPENDIX D 

Orthogonality 0/ the Modes 

Let q, p and q\, p be any two solutions of the paraxial ray equations, 
equations (11a and b), and assume that the matrix PQ~\ where 

Qm[qQ, (125a) 

P^[pp], (125b) 

is symmetric. 

An infinite set of solutions of the parabolic wave equation has been 
obtained in the main text in the form [equation (69)] 

+ mn (r, z,Q,P)= \Q I" 1 exp (-i | fPQ" v)tf mn (x ; r), (126) 

where H mn denote the Hermite polynomial in two variables xi , X2 

X - Q*~ l r, (127) 

associated with the quadratic form x>Xi where 
v = jkQ*(PQ-> - P*Q*-')Q* 
= JQ-'Q*, (128) 

J - - J** ~ Pq * Y ~ ~™*\ (129) 

_qp* - PQ* ~qp* - pq*_ 

Let us now impose on the rays the additional condition 

(<R; (R*) = qp* - p§* = (130) 

and assume that the diagonal terms of J are positive. Since, as pointed 
out in the main text, the two rays need be denned only to within con- 
stants, they can be normalized in such a way that J is the unit matrix. 
In that case we have, from equations (127) and (128) 

(131) 



(132) 



"X = x 
Consequently 

H*. n .(x; ») = H m . n .( x *; "*) = G m . n .( X ; v) (133) 

where we have introduced the adjoint polynomials G mn defined by 

G nn { X ;v) ^ H mn {v X ;v- 1 ). (134) 
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The orthogonality condition for two solutions \J/ mn and ^ m -„< [equation 
(126)] can now be written in the form 

= I Q~ l Q* | 5 //*" exp (-tevx)H mn (x',v)G m , n ,( X ','>) d\, 

= 2irm\n\ if m' = m and n' = n, 

= if m' ?^ m or n' ?± n. (135) 

The biorthogonality property 32 of the polynomials H mn and G mn has 
been used in equation (135). 
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